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Abstract 

We reconsider the high energy resummation of photoproduction, electroproduction and 
hadroproduction cross-sections, in the light of recent progress in the resummation of per- 
turbative parton evolution to NLO in logarithms of Q 2 and x. We show in particular that 
the when the coupling runs the dramatic enhancements seen at fixed coupling, due to in- 
frared singularities in the partonic cross-sections, are substantially reduced, to the extent 
that they are largely accounted for by the usual NLO and NNLO perturbative corrections. 
This leads to a novel explanation of the large X-factors commonly found in perturbative 
calculations of hadroproduction cross-sections. We give numerical estimates of high energy 
resummation effects for inclusive S-production, inclusive jets, Drell-Yan and vector boson 
production, along with their rapidity distributions. We find that resummation modifies 
the S-production cross-section at the LHC by at most 15%, but that the enhancement of 
gluonic VF-production may be as large as 50% at large rapidities. 
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1. Introduction 

At the LHC we hope to separate with confidence a tiny fraction of interesting events 
from an overwhelming background of collisions involving gluons carrying only a small 
fraction of the momentum in the beams. The success of this enterprise depends crucially 
on our ability to control high energy logarithms in perturbative QCD at leading twist when 
we calculate inclusive cross-sections. Currently however no reliable calculations of high 
energy resummation corrections to any hadronic process have been made. The purpose of 
this paper is to remove the one remaining obstacle to performing such calculations. 

In recent years there has been considerable progress by several groups in understanding 
the resummation of parton evolution, so that we now know how to simultaneously resum 
all collinear and small- £ logarithms at NLO. This programme depended on several key 
ingredients: /c^-factorisation [1,2], NLLx corrections [3,4,5], the recognition of the need to 
simultaneously resum collinear, anti-collinear and high energy logarithms [6,7], the use of 
high energy duality to achieve this [8,9], and the understanding of running coupling effects 
[10,11,12,13]. It is now possible to perform precise and reliable calculations of small- £ 
resummation corrections to parton distribution functions. 

Hadroproduction processes have received much less attention [14,15,16,17]. The reason 
for this is partly their additional kinematic complexity, but also because of a difficult 
conceptual problem standing in the way of reliable results. This problem relates to the 
infrared singularity which appears when two gluons collide at high energy, due to the 
possibility of all the energy going into a timelike gluon which may then go almost on-shell 
[15]. Though it has been understood for some time that this singularity might produce 
substantial enhancements of hadronic cross-sections, it has been difficult to make reliable 
predictions, particularly when the coupling runs [16]. This is the problem we resolve in this 
paper. We will find that the singularity is less dangerous than naive arguments suggest, 
that the enhancements it produces at high energy are modest, and in fact may be well 
approximated by the NLO and NNLO perturbative results. 

The structure of the rest of the paper is as follows: in chapter 2 we summarise the 
main ideas used in the resummation of high energy logarithms in order to set the scene 
and fix notation, and then explain the difficulties encountered in applying these ideas to 
the resummation of hadronic cross-sections due to infrared singularities. In chapter 3 we 
consider the simpler scenario of photoproduction and electroproduction processes, and in 
particular how there infrared singularities may be dealt with using an exponentiation trick. 
We apply this trick to the inclusive photoproduction of 66-pairs, providing quantitative es- 
timates of resummation effects. We then in chapter 4 move on to the more interesting 
case of hadroproduction, construct the gluon-gluon luminosity, describe the singularity 
structure of the partonic cross-sections, and show how the same trick used in photopro- 
duction works here also. We provide generic estimates for resummation effects in various 
hadroproduction processes at the Tevatron, LHC and a notional VLHC, and consider in 
detail the particular case of hadroproduction of bb pairs. We also consider the stability 
of the resummed perturbative expansion. In chapter 5 we consider how we may compute 
rapidity distributions in this framework, and offer estimates of resummation corrections 
for the rapidity distributions of bb pairs and W bosons at LHC. Finally in chapter 6 we 
summarise our results, and suggest several directions for future work. 
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2. High Energy Singularities 

2.1. High Energy Factorization 

We consider electron-hadron, photon-hadron and hadron-hadron processes in which 
Q is the hard (transverse) scale (for example the photon virtuality, a heavy quark mass, 
or the invariant mass of some particular particles in the final state), S the square of the 
centre-of-mass energy, and p = Q 2 /S. The dimensionless cross-section E = Q 2 a is a 
function of p and Q, the scale of Q being set by Aqcd- For hard processes Q 3> Kqcd, 
while for high energy processes S ^> Q 2 , so p <C 1. 

In an electroproduction or photoproduction process, if x and k are the longitudinal 
and transverse momenta of the struck parton, the square of the centre-of-mass energy of the 
hard process is s = xS, and the dimensionless hard cross-section E 7J - = Q 2 a 1 j is a function 
of Q 2 / s = p/x, k/Q and p/Q, where p is the factorization and renormalization scale (here 
set equal) , and j labels the struck parton. The (unintegrated) parton distribution function 
Tj depends only on i, k 2 and p 2 . Factorization (or more specifically "/cT-factorization" 
[1,2,14,15]) is then the statement that the photon-hadron cross-section may be written as 

S 7h (p,Q)= £ f^r^^^a.ip 2 ))^^, 2 ), (2.1) 

j=g,q,q p 

up to terms which vanish as inverse powers of the hard scale Q. 

For a purely hadronic process, the centre-of-mass energy of the hard process is s = 
X1X2S, where x\ and x^ are the longitudinal momentum fractions of the colliding partons 
ji and J2 in hadrons hi and h%. The dimensionless hard cross-section E Jlj2 = Q 2 o~j 1 j 2 is 
then a function of Q 2 /s = p/x\X2, ki/Q, k 2 /<5 and p/Q, so factorization is the statement 
that the hadron-hadron cross-section may be written as 



v ( n s ^ f 1 dx 1 [ 1 dx 2 fd 2 ^ fd 2 k 2 



3i,32=g,q,q 



P *i J P ^2 J «m J 7rk| ^ ^ 

again up to terms which vanish as inverse powers of the hard scale Q. The dependence on 
p will be suppressed in what follows: in practice we will take p = Q. 

2.2. Gluon dominance at high energy 

The range of x and Q relevant at lepton-hadron and hadron-hadron colliders (specif- 
ically HERA with VS = 320GeV, the Tevatron with VS = 1.8TeV, the LHC with 

= 14TeV and a notional VLHC with y/S = 200TeV) is shown in fig.l. It is obvi- 
ous from the figure that over most of the kinematic reach of these machines x is small, and 
thus small x logarithms are potentially large: only for processes at the highest scales (for 
which the cross-sections are correspondingly small) can the small x region be altogether 
excluded. 
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Figure 1: the kinematic plane for electro-hadron and hadron-hadron collisions: x is the 
fraction of longitudinal momentum in the parton, and Q is the hard scale in GeV (which 
for photoproduction and hadroproduction is the invariant mass of the particles produced 
in the hard process). Shown are the regions accessible in fixed target experiments 
(orange), at HERA (red), at the Tevatron (green), at LHC (blue) and at a notional 
VLHC with S = 200TeV (grey). The dashed lines for the hadron colliders show the 
central rapidity region, where x\ = xi = x: when the rapidity is nonzero, the values of 
x\ and X2 may be read off for a given Q by choosing points symmetrically placed about 
this line. The solid black diagonal line shows the points at which L z = L' z , and thus 
the logarithms of Q 2 are as important as the logarithms of x (see eqn.(6.1) below): the 
dotted black line is determined by L z = 0.1L' Z , while the heavy black line is determined 
by L z = 10L' Z . 
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(a) (b) (c) 



Figure 2: the quark induced Drell-Yan process, or vector boson production (a) the 0(a) 
qq annihilation process and (b),(c) the 0(aa B ) qg (or qg) fusion process, with initial 
and final state radiation respectively. 




(a) (b) (c) 



Figure 3: the gluon induced Drell-Yan process, or vector boson production: the three 
main classes of contribution (a) with initial state radiation, (b) with both initial and 
final state radiation and (c) with final state radiation only. All these processes are 
formally 0(ota 2 s ). 

High energy processes thus usually involve collisions of small x partons, and at small 
x and high Q 2 it is now well known that all singlet parton distributions show a steep rise, 
that of the gluon distribution G(x, Q 2 ) being driven by the nonabelian splitting g — > gg, 
and that of the singlet quark distribution Q(x, Q 2 ) by g — > qq [18,19]. Because the latter 
process is 0(a s ), at small x and high Q 2 the singlet quark distribution is always smaller 
by a power of a s (Q 2 ) than the gluon, i.e. Q(x,Q 2 ) ~ a s (Q 2 )G(x, Q 2 ): at small x most 
partons are gluons. 

This means that when we compute high energy partonic cross-sections we can no 
longer rely on the simple counting of powers of a s as we do when Q and S are comparable. 
Consider for example a Drell-Yan type of process: formally the LO process is qq annihila- 
tion, fig. 2(a), and is 0(a), while the qg (and qg) scattering processes fig. 2(b) and fig. 2(c) 
are both NLO, i.e. 0(aa s ), and the gg processes fig. 3 are all 0(aa 2 ) and thus NNLO. 
However at high energy, if we take into account the relative suppression of the quark rel- 
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ative to the gluon, all of these contributions are in practice 0(aa 2 ), and thus should be 
considered as leading order. Subleading contributions to the qq and qg contributions in 
fig. 2 are known [20]: the subleading contributions to the gg processes fig. 3 are not known 
at present, since formally they would be NNNLO. 

In this paper we shall only consider processes with gluons in the initial state. In 
practice this means that we simply drop the summation over partons in eqns.(2.1),(2.2). 
This simplification is sufficient to discuss most of the issues in high energy resummation, 
but will of necessity mean that our numerical results are estimates rather than calculations 
of complete cross-sections. 

2.3. Double Mellin Transforms 

The convolutions over x and k in the factorizations (2.1) and (2.2) may be undone by 
taking Mellin transforms with respect to p and Q. Explicitly if we define 

S 7h (iV, M) = £% N f^f^) " M S 7 h(p, Q% (2.3) 
while for the hard cross-section 

CW M) = /fV/§(|) M E w (4), ,2.4) 

then if the double Mellin transform of the unintegrated gluon distribution Q = T g is 

G(N,M)= —x N -tf -ro Q(x,k 2 ), (2.5) 



k 2 \A 2 J 

the factorization (2.1) becomes simply algebraic: the purely gluonic contribution to the 
transformed cross-section is 

E 7h (iV, M) = C(N, M)Q(N, M). (2.6) 

It is easier to make contact with phenomenology if instead of working with the unin- 
tegrated gluon distribution we define the integrated distribution 

r Q2 rlk 2 

G(x,Q 2 ) = j^ -^Q{x,k 2 ). (2.7) 

This is the distribution that would be proportional to the physical cross-section if the hard 
process were pointlike, i.e. the hard cross-section was simply proportional to <S)(Q 2 — k 2 ). 
Then in Mellin space 

G(N,M) = M~ 1 g(N,M), (2.8) 

and eqn.(2.6) becomes 

£ 7h (iV,M) = C(N,M)G(N,M). (2.9) 
where C(N, M) = MC(N, M). 
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To recover the physical cross-section we invert the two Mellin transforms (2.3): 

/ io ° rIN f io ° rIM 

7^/ N J i |^ M C(iV,M)G(iV,M), (2.10) 

where for clarity in future discussions we have defined 

e = logl/p, t = logQ 2 /A 2 . (2.H) 

Here N and M are both complex variables: the contours in the integrations over N and 
M keep just to the right of the singularities near N = and M = 0. The contour in N 
is always closed on the left (since £ > 0): the contour in M is closed on the left in the 
ultraviolet (t > 0), but on the right in the infrared (t < 0). 
For hadronic processes we proceed similarly: defining 

-HIM M M\ K M'ki/kfxMi /-d 2 k 2 /k|xM 3 k, kj 

WM " Mi)= /,//^M / s »("-q'o ) - (2 ' 12) 

the factorization formula (2.2) becomes 

£ hh (iV, M u M 2 ) = H(N, M u M 2 )G{N, M 1 )G(N, M 2 ). (2.13) 
where H(N, Mi, M 2 ) = M 1 M 2 H(N, Mi, M 2 ). The hadronic cross-section is thus 

£ hh (p,Q)= r^e* N P^^e^ + ^^(iV,Mi,M 2 )G(iV,Mi)G(iV,M 2 ), 

(2-14) 

and we have to deal with functions of the three complex variables N, M\ and M 2 integrated 
along three contours. Clearly H(N, Mi, M 2 ) = H(N, M 2 , Mi). 

2.4- Duality at fixed coupling 

Perturbative expansions of the hard cross sections are contaminated by logarithms of Q 
and p, corresponding in Mellin space to inverse powers of M and iV respectively. This may 
be seen directly by considering Laurent expansions around M = iV = of the integrands 
of (2.10) and (2.14): every extra inverse power of M or N must be compensated by a 
positive power from expansion of the exponential, yielding an extra factor of t or £. These 
logarithms are at most single logarithms, in the sense that in perturbation theory there is 
at most one extra logarithm of each type whenever there is an extra power of a s : a typical 
contribution to the integrand is a l s M~ m N~ n where m,n < I. To obtain meaningful results 
in perturbation theory these logarithms must be resummed (at LO, where m + n = 2/, 
NLO, where m + n = 21 + 1, etc.) and factored into the gluon distribution G(M, N). In 
particular without resummation the GLAP splitting function is unstable in the small x 
region, while the BFKL kernel is unstable in both the collinear and anticollinear regions, 
which means that for reasonable values of a s it is also unstable in the small x region [6] . 
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Resummation of the transverse and longitudinal logarithms t = log<5 2 /A 2 and £ = 
logl/p proceeds by solution of the GLAP and BFKL equations respectively: 

f = /'! P (^(Q 2 ))G(*,Q% (2.15) 
where P is the gluon splitting function, and 



dg^ r 



^/C(^,a s (fc 2 ))£(z,fc 2 ), (2.16) 



where /C is the (angular averaged) BFKL kernel. Taking double Mellin transforms both 
equations simplify to algebraic equations, since the convolutions reduce to products and 
the derivatives to factors of M and N respectively: 

MG(N, M) = G (N) + 7 (iV; a s ) G(N, M), 

(2.17) 

NG(N, M) = Gq(M) + M x(M, a s )M G(N, M), 
where 7 and x are the Mellin transforms of the respective kernels 

™x N P(x;a 8 ), X (M;a s ) = ^ ^ (^) /C(^, a s (fc 2 )), (2.18) 

and Gq(N), Gq(M) are (nonperturbative) boundary conditions. The transform of the 
BFKL equation has been written in terms of the integrated distribution using eqn.(2.8). 

The main complication in the derivation of eqns.(2.17) from eqns(2.15),(2.16) is the 
running of the coupling: in Mellin space a s (t) becomes an operator 

= ( 2 - 19 ) 

This leads to various difficulties, discussed at length in refs[ll,12,21,22,23]. In the remain- 
der of this section we sidestep this issue, and consider the coupling to be fixed. 

At fixed coupling, a s — > a s , and the evolution equations have the simple algebraic 
solution 

G(N, M) = — L^Go(iV) = 1 -—-Go{M). (2.20) 

M-7(iV;a s ) N-x{M;a s ) 

Since the leading twist perturbative singularities can always be factorised from the singu- 
larities in the nonperturbative boundary conditions, the poles in the perturbative factors 
must coincide: at the pole 

M = 7 (iV;a s ), N = X {M-a s ). (2.21) 
The functions 7 and x must thus satisfy the consistency conditions 

M = 7 ( X (M; a a ); a 8 ), N = xh(N; a s );a s ), (2.22) 
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i.e. as functions of their main arguments 7 = x -1 , x = 7 _1 - These are the duality relations 
[8,14,24,25]: x determines the high energy (N = 0) singularities of 7, just as 7 determines 
the collinear (M = 0) singularities of x- F° r example if 7 ~ ot 8 /N, x ~ ol s /M and vice 
versa. Conservation of longitudinal momentum implies that 7(1, a s ) = to all orders in 
perturbation theory: duality then tells us that x(0; a s ) = 1 to all orders in perturbation 
theory, i.e. that the collinear (M = 0) singularities in the expansion of x(M; a s ) in powers 
of a s resum to unity [9]. 

Using these ideas, and the symmetry under the exchange of Q and k of the BFKL 
kernel (which in Mellin space translates into a symmetry under M — > 1 — M [7]), it is 
possible to use the the known results for the GLAP anomalous dimension 7(iV; a s ) at LO 
and NLO in powers of a s to resum the collinear (Q 2 3> k 2 , thus M = 0) and anticollinear 
(Q 2 < k 2 , thus M = 1) singularities in the BFKL kernel x(Af ; a s ) at LO and NLO. Indeed 
the resummation of these singularities is essential to obtain a meaningful expansion of the 
kernel for reasonable values of a s . The resummed kernel then in turn through duality gives 
an anomalous dimension in which the high energy N = singularities are also resummed. 
This resummed anomalous dimension can then be used to evolve an initial (integrated) 
gluon distribution at small x [9,26]. 

The small-x behaviour of the fixed coupling anomalous dimension (or rather its as- 
sociated splitting function) is given by the behaviour around the minimum of %(M, a s ): 
fixed coupling duality implies that this leads to a square root branch cut at iV = c(a s ) = 
x{\]ot s ) : at which the anomalous dimension rises to one half: 

1 \N — c(a s ) 

V i K ( a °) 

where K,(a s ) = x"(^, a s ) is the curvature at the minimum. This cut in turn gives rise to 
the famous x~ c ^ asS) growth in the splitting function. 

Fixed coupling duality may also be used to resum the high energy logarithms in hard 
cross-sections. Since all the collinear and high energy logarithms have now been absorbed 
into the integrated gluon distribution, the hard cross-sections are regular in both iV and 
M close to the origin, and may thus be Taylor expanded: for example 

00 00 

C(iV,M;a s ) = E cl(N)a l s M m , 

1=1 m=0 

H(N,M 1 ,M 2 ) = f^ f; £ h l mim2 (N)a l s MrM2 

1=2 mi =0 m2=0 

where c l m (N) and h mim2 (N) = h m2Tni (N) are also regular in the neighbourhood of iV = 0. 
The only singularity close to the origin is thus that in G(N, M), i.e. (2.20), and this may 
be used to perform one of the photoproduction inverse Mellin transforms (2.10), or two of 
the hadroproduction inverse Mellin transforms (2.14). The usual procedure is to perform 
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(2.23) 



(2.24) 



the integrals over M in this way, leaving the single integral over TV to be done numerically: 



/lOO J AT 
« ^)C(JV, 7 (iV; a.); « s )G (iV), 
-too lm 



-too 
cioo 



e^V^^tf (AT, 7 (iV; a s ), 7 (AT; a s ); « s )G (iV) 2 . 

-too 2ya 



(2.25) 



Note that since 7 (iV; a s ) contains resummed poles in N, (i.e. terms of the form 
a l s N~ n , with n < I) this procedure effectively resums the high energy singularities in the 
collinearly factorised hard cross-sections, obtained by expanding C(N, 7 (iV; a s ); a s ) and in 
H(N, 7 (iV; a s ), r y(N; a s );a s ) respectively in powers of a s [2,15,25]. Conversely, comparing 
eqns.(2.25) with the more conventional expressions 



/tOO 7 AT 

^ N e^ N ^c(N;a 8 )G (N), 
-too lm 

/ io ° dN 
^l e ^e 2 ^ N ^f(N;a s )Go(Nf, 
-ioo 2nl 

obtained in collinear factorization, we see that 

oo 

al s c i( N ) = C(N, ^(N; a s ); a s ), 

1=0 

oo 

f(N; a a ) = al sfi( N ) = H ( N i 7(^5 7^; a s ); a s ). 



(2.26) 



(2.27) 



1=0 



Thus by substituting the fixed order expansion 7(iV; a s ) = Xl^o a sli(N) into the expan- 
sions (2.24) of the right hand side of (2.27) the high energy singularities of the fixed order 
hard cross-sections ci(N) and fi(N) may be economically computed. This makes matching 
to the fixed order calculations particularly straightforward. 

This is the procedure used in most previous studies [2,15,16,17,25,26] of the effects 
of high energy resummation on high energy cross-sections. It leads to a strong growth in 
cross-sections, particularly hadronic cross-sections, due to infrared singularities at positive 
values of M which, unlike the collinear singularities near M = 0, have not been resummed 
into the perturbative evolution. It is the nature of these infrared singularities, and the 
correct treatment of them, that is the main subject of this paper. First however we must 
consider how this simple picture is modified when the coupling runs. 

2.5. Duality at running coupling 

When the coupling runs the evolution equations (2.17) become differential equations 
in M, since the running coupling transforms to a differential operator (2.19). This operator 
commutes with iV but not with M: 

[a s ,M} = p(a 8 ), (2.28) 
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where (3(a 8 ) = — ol 2 s (3q-\ is the QCD beta-function. This means that the BFKL operator 

x(M; a s ) must be defined very carefully: different orderings of a s and M will give different 
results, because the arguments of the couplings will be different [22]. 

It was show some time ago both by saddle point arguments for the Mellin inversion [27] 
and by explicit solution of the M-space differential equation [28] that when the coupling 
runs, the naive duality (2.22) is modified by terms proportional to (3: for example at NLLx, 
writing X (M, a s ) = a sXo (M) + a 2 sXl (M) + ■■■, 

Xi^Xi + A)^, (2.29) 

where the primes denote derivatives with respect to M. In ref.[22] it was shown that the 
fixed coupling duality (2.22) may be extended formally to relations between the operators 
7(-/V;ai s ) and x(M;a s ), which may then be used to generate systematically running cou- 
pling corrections using a purely algebraic algorithm, order by order in perturbation theory. 
This technique was used recently [23] to compute an estimate of the leading twist BFKL 
kernel X2 at NNLLx. 

When the coupling runs it becomes important to specify carefully the factorization 
scheme, since a redefinition of the gluon distribution G{N, M) — > Z(M)G(N, M) changes 
x(M, a s ) by a commutator: 

X (M, a s ) - X (M, a 8 ) + Z~\M)\ X {M- a s ), Z{M)\. (2.30) 

It is thus possible to shuffle the running coupling terms between the evolution and the 
hard cross-section. In this paper we will use the QoMS scheme [29,30,31], a variant on the 
MS scheme in which all the running coupling terms are absorbed into the evolution. The 
two schemes are equivalent in fixed order perturbation theory at LO, NLO and NNLO, 
but begin to differ at NNNLO. 

Note that provided Z(M) is regular at M = the second term on the right hand side 
of (2.30) is subleading compared to the first, since the commutator necessarily introduces 
an extra power of a s (see (2.28), and note that (5{a s ) = —floa 2 . + 0(a^)). This means that 
if we work at NLLx in the resummation of the singularities, it is only necessary to consider 
the hard cross-sections C(N, M; a s ) and H(N, Mi, a s ) at LLx: terms at NLLx in the 
cross-section are NNLLx in the evolution [2,28]. Similarly to this order we may ignore the 
running of a s in the hard cross-section, setting a s — > a s up to subleading terms. 

Away from the region of very small x the running coupling corrections to naive duality 
are small. However in the small x limit, they become very large. This may be seen 
immediately from (2.29): at small x and fixed coupling, the integrals over M are dominated 
by the minimum of x(M, ct s ), and thus by M ~ |, ie the region in which xb(M) ~ and 
the correction (2.29) becomes large. It was thought at one time that these terms become 
so large that they drive instabilities in the gluon distribution leading to negative cross- 
sections [10,30,32,33]. However these instabilities are due to diffusion into the infrared: if 
the singularities are resummed, they may be factorised into the (nonperturbative) initial 
distribution, resulting in stable evolution [11,12,13]. 

This further resummation of the running coupling singularities at M = | is accom- 
plished by expanding around the minimum of x, solving for G(N, M) and then performing 
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1 10 100 1000 10000 



Q/GeV 

Figure 4: The gluon distribution G(x, Q) plotted against Q in GeV for x = 
1CT 10 , 1(T 9 , ... ,0.01,0.1 (from top to bottom). The blue curves are evolved with the 

NLO resummation described in the text: the black and red curves are with LO and 
NLO GLAP evolution respectively. 
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the inverse Mellin with respect to M exactly, rather than using a saddle point expansion. 
The simplest version of this argument gives rise to an Airy function [21]: a full calcula- 
tion [12] summing up the leading singularities in (3q gives the running coupling resummed 
anomalous dimension 

7s (A,a s ) = - 2 - (3 a s + jK2BmasA) - iy (2-31) 

where A and B are simple functions of a 8 and N computed from the value of x(M, a s ) 
and its curvature at the minimum, and K v {x) is the Bateman function. The small x 
behaviour resulting from (2.31) is qualitatively different from that obtained with the fixed 
coupling result (2.23), since the cut is replaced by a simple pole located at N = cb(cc s ), 
with cb{ol s ) given by the rightmost zero of the Bateman function. Since cb(cx s ) is rather 
less than c(a s ), the x~°^ as > growth of the splitting function at small x is softened by the 
running coupling effects to x~° B ^ as > [12]. However since the new singularity is now a pole 
and not a cut M can now grow indefinitely rather than saturating at M = |: in effect the 
region between M = and M = \ is stretched to infinity so that the effective %(M; a s ) 
(the naive dual of 7(A; a s )) is analytic for all M > 0. 

By combining the running coupling resummation with a collinear and anticollinear 
resummation of x(M; a s ) using running duality, a completely resummed anomalous di- 
mension may be computed, and used to resum both high energy and collinear singularities 
in the gluon distribution at LO and NLO [12,13]. Such a NLO resummed gluon distribu- 
tion is plotted in fig. 4 for x down to 10~ 10 and Q up to 50 TeV. The initial distribution 
was chosen to be proportional to x (1 — x) 5 at 2 GeV, with the first moment (i.e. the 
integrated longitudinal momentum) normalised to unity. We take a s (mz) — 0.118, a 8 (Q) 
evaluated using standard two loop running with thresholds at mj and m t . For compari- 
son we also plot the same distribution evolved using the LO and NLO GLAP anomalous 
dimensions. All evolution is performed with n/ set to zero in the anomalous dimensions, 
consistent with our suppression of all quark induced processes. It can be seen from the 
plot that the effect of the resummation is modest, even over such a wide kinematic range, 
and is generally such as to smoothly soften the growth at small x and large Q 2 . 

2.6. Soft Singularities in Hard Cross-sections 

All that remains to be done is to combine this resummed gluon distribution with a 
hard cross-section. However when the coupling runs it is no longer so easy to justify the 
use of the pole approximation (2.20) to perform the integrals over M in (2.10) and (2.14) 
as we did at fixed coupling, eqn.(2.26), since for the gluon distribution itself we know 
that at small x the pole approximation no longer gives the correct asymptotic behaviour. 
Moreover, since when the coupling runs the resummed anomalous dimension 7 (iV, a s ) rises 
indefinitely as N decreases, the position M = 7(77,, a 8 ) of the pole in the evolved gluon will 
also grow large, and thus eventually the contour (which is always to the right of this pole) 
will become entangled with the various infrared singularities in the hard cross-sections. 

Consider first the photoproduction or electroproduction hard cross-section C(N,M), 
which may be computed at leading order from the graph shown in fig. 5. The hard scale Q is 
then the mass of the quark (for photoproduction of heavy quarks), or the photon virtuality 
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Figure 5: Photoproduction of heavy quarks. The photon is on-shell, while the gluon 
is off-shell by an amount k 2 , in order to calculate the impact factor C(N,M). In 
electroproduction the photon would also be taken off-shell. 



(for electroproduction). The incoming gluon has virtuality k 2 . At very high energies 
S ^> Q 2 , relevant values of iV should be small, so can in the first instance be neglected. 
Then considered as a function of M, C(0,M) is regular near M = (since the collinear 
singularities, with Q 2 3> /c 2 , have already been absorbed in the gluon distribution), but has 
poles at M = —1, —2, . . . and M = 1, 2, 3, . . .. The poles at negative values of M correspond 
to higher twist singularities: they lead to power corrections in inverse powers of Q 2 , which 
are not relevant here. The poles at positive values of M correspond instead to process 
dependent infrared (anticollinear) singularities with Q 2 <C k 2 . It is these singularities that 
enhance the cross-section at high energy. 

When the coupling is fixed, the anomalous dimension for the gluon evolution saturates 
at the cut, and thus at high energy the dominant contribution to the M integral is from 
the region M = |. The resummation K- factor may be then be estimated at high energy 

to be C ^qq^ = ngf- ~ 4.8, a significant enhancement [15]. However when the coupling 
runs, there is no such saturation, and the pole approximation to the integration over M 
breaks down at high energy, since the pole at M = 7(iV, a s ) in the evolved gluon can 
then approach the pole at M = 1 in the hard cross-section, and the contour gets pinched 
between them. We must then attempt to perform the integration exactly, and afterwards 
factor out the gluon distribution: resummation of the hard cross-section is thus rather 
more difficult than it was at fixed coupling. 

This difficulty becomes much worse when we consider hadroproduction. Once again 
we focus on the singularities of the impact factor H(0, M%, M2) in the double-Mellin plane 
of Mi and M2. This may be obtained from the graphs in fig. 6 (hadroproduction of heavy 
quarks, so Q is the quark mass), or fig. 6 and fig. 7 (inclusive jets, so Q is the transverse 
momentum of the jets). Again the region around the origin is regular, and there are 
higher twist poles at Mi, M2 = — 1, —2, . . . and infrared singularities at Mi, M2 = 1, 2, . . .. 
However now we also have lines of infrared singularities at Mi +M2 = 1, 2, . . .: for example 
the heavy quark hadroproduction impact factor 

g(M 1 ,M 2 )~ a ^ (1 _ Mi 1 _ M2)3 , (2.32) 
when Mi + M2 ~ 1 [15,16,27]. It is easy to see why this occurs by noting the form of 
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(a) (b) 



Figure 6: Hadroproduction of heavy quarks: (a) the "abelian" diagram, with initial state 
radiation and (b) the "nonabelian" diagram with final state radiation. Both incoming 
gluons are taken off-shell by amounts k\ and k% in order to calculate the impact factor 
H(0, Mi, M2). The same graphs give hadroproduction of quark jets. 




(a) (b) 



Figure 7: Hadroproduction of gluonic jets: (a) the "abelian" diagram with initial state 
radiation and (b) the "nonabelian" diagram with final state radiation. Both incoming 
gluons are taken off-shell. The two final state gluons are on-shell, with a specified pr- 
The four gluon vertex diagram is not shown, since it is nonsingular. 

the Mellin transform eqn.(2.12): an infrared singularity when Q 2 -C k 2 ,k 2 will become a 
singularity at My + M 2 = 1, 2, . . .. 

The degree of this singularity will depend on the nature of the infrared singularities of 
the individual graphs making up the process. The graphs with initial state radiation, fig. 6a 
and fig. 7a, both have a single (anti)collinear singularity similar to that in photoproduction 
fig. 5. However the singularity structure of the nonabelian graphs with final state radiation 
fig. 6b and fig. 7b is more complicated. Firstly there is a soft singularity when the internal 
(timelike) gluon goes on-shell: the denominator of the propagator is s = X1X2S— (ki+k^) 2 , 
and the singularity arises when X1X2S ~ (ki + k 2 ) 2 <C k 2 ~ k| [15]. On top of this there 
is the usual collinear (heavy quark production) or soft (jet production) singularity for the 
emission of the final state partons, so altogether we find the triple pole eqn.(2.32). Note 
that this kind of singularity is thus generic to most gluon-gluon hard cross-sections. In 
fact it is the dominant singularity, providing most of the cross-section at very high energy. 
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To evaluate the cross-section, we must perform the double Mellin inversion eqn.(2.14). 
When the coupling is fixed, M\ and M 2 both try to saturate at M\ = M 2 = |, and 
thus touch the singularity. This produces a strong enhancement of the cross-section [15], 
stronger even than the BFKL growth. If M\ and M 2 can grow beyond one half, as one 
expects at running coupling, the contours are pinched, and the cross-section seems to grow 
as fast as p~ 2aa [16]. However this dramatic growth is also very unstable, in particular to 
N dependent corrections, since these split the triple pole into a separate double pole (at 
Mi + M 2 = 1 + N) and a single pole (at Mi + Mi = 1). This instability is a sure sign that 
the resummation is not under control. Indeed, when the coupling runs, the whole line of 
singularity become accessible (since both Mi and M 2 may become large), and furthermore 
becomes entangled with the collinear regions (since at Mi ~ 0, M 2 ~ 1 and vice versa). 

So somehow we need a more reliable way of computing the integrals over Mi and M 2 . 
Since the singularity (2.32) does not factorise into a function of Mi times a function of 
M 2 , we actually need to perform both integrals simultaneously, and then factorise out the 
gluon distributions afterwards to obtain the resummation of the hard cross-section. 

In the rest of this paper we will show how a simple trick may be used to solve this prob- 
lem, first for photoproduction and electroproduction (Section. 3) and then for hadropro- 
duction (Section. 4). We consider processes with a single incoming gluon first because they 
are simpler: the more important hadroproduction processes will then be dealt with by a 
straightforward extension of essentially the same idea. 



3. Photoproduction and Electroproduction 

3.1. Evaluating the Cross-section 

Consider first a photoproduction or electroproduction cross-section at high energy, and 
thus small p, but still far from the high energy limit. Then we might expect M to remain 
"small" in some sense, and consider approximating the hard cross-section C(N, M; a s ) by 
the first few terms in its Taylor expansion (2.24). The cross-section (2.10) may then be 
written 

/too jnr 
e^ N E lh (N,t), (3.1) 
-ioo 2nl 



where 



E 7h (iV, t) = a s (t) cl(N) / 7j — 7 e tM M m G(N , M), (3.2) 



where we have kept only the leading term in the expansion in a s , and (optimistically) 
changed the order of integration and summation over m, in order to do the integrals term 
by term. However for the first term this is an integral we already know: 

G(N,t)= I |-e tM G(iV,M), (3.3) 



— too 



and consequently for all m = 0, 1, 2, 



100 rJM 8 m 

—e tM M m G(N,M) = —G(N,t). (3.4) 
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1000 10000 

Q/GeV 

Figure 8: The first (left) and second (right) derivatives with respect to t of the gluon 
distributions G(x, Q) plotted in fig. 4, normalised to G, plotted against Q in GeV, for 
x = 1CP 10 , 10~ 9 , . . . 0.01 (from top to bottom). Again the blue curves are with NLO 
resummed evolution, the red with NLO GLAP evolution. 



However the general solution of the evolution equation (2.15) is simply 

G(N,t) =exp(^ dt' 7 (N,a s (t'))^Go(N), (3.5) 

so all the partial derivatives in (3.4) may be evaluated in terms of the anomalous dimension 
and its derivatives: 

^G(N,t)= 1 (N,t)G(N,t), 

^G(N,t) = ( 7 2 + j)G(N,t), (3.6) 
d 3 

W G(N, t) = ( 7 3 + 3 7 7 + l)G{N, t), . . . 
where the dot denotes partial derivatives with respect to t. 



16 



Note that 7 is formally subleading compared to 7, since 



s 



so every time we differentiate with respect to t we add a power of a s . So up to subleading 
terms partial derivatives of G(N, t) with respect to t simply result in powers of 7(iV; a s (t)), 
and combining eqn.(3.2) with (3.4) we find 

00 0777, 

E 7h (iV,t) = a s (t) cl m(N)-^G(N,t) = a s (t)C(N, 7 (iV, a s (t)))G(N, t). (3.8) 

m=0 

So provided we are not close to a pole, we find the same result when the coupling runs as 
we found at fixed coupling case through the pole dominance argument (2.26). 

Formally this method works provided the series (3.8) converges. To see whether this is 
likely in practice, we plot in fig. 8 the first and second derivatives of the gluon distributions 
fig. 4, normalised to G: indeed the derivatives are much less than one for all but the lowest 
values of x and Q. Moreover the second derivatives are much smaller than the first, 
suggesting that the series converges rather rapidly. In addition derivatives for the NLO 
resummed gluons are rather smaller at low x and Q than those for the NLO GLAP gluons, 
essentially because they evolve rather more slowly. It follows that at all but the smallest 
x and Q we may use Taylor expansions in M, and we will show repeated examples of this 
in what follows. It is worth noting in parenthesis that Taylor expansion in M is little 
different at high energy to the usual fixed order expansion of the hard cross-section (2.27), 
since it is only through the M dependence of C(N, M; a s ) that small N singularities can 
be introduced, and at high energy it is these that dominate the cross-section. 

3.2. Resumming an Infrared Singularity 

Now consider what happens when C(N, M; a s ) has a pole near M = 1, so that the 
Taylor expansion around M = has radius of convergence one. This will generally be the 
case: such a singularity corresponds to logarithms of k 2 /Q 2 as k 2 — > 0. We will then need 
to do integrals of the form 

/ io ° rlM 1 
^ilT^W^W), (3.9) 

for n = 1, 2, . . .. We can perform these integrals exactly using a simple trick [34]: we write 

_ i_ - = I r dTT n - x e-^- M \ (3.10) 
(1-M)» nlj K ' 

thereby exponentiating the dependence on M and transferring possible singularities in the 
complex variable M to singularities in the integration over the real variable r. Substituting 
(3.10) into (3.9) the M integration is now indeed trivial: 



Zj 7 h 



/ lo ° rlM 1 f°° 
^L-l drr n - l e- T+M{t+ ^G(N,M) 
-too 2m n] J0 (3 n ) 

dr r n ~ 1 e~ T G(N, t + r), 
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where in the second line we exchanged the order of the integrations over M and r and 
performed the integration over M. Since when the coupling runs the growth of G(N,t) 
with t is rather gentle (no greater than a power of t), the r integral converges for all t. At 
small x the M = 1 singularity thus always leads to a modest enhancement of G(x, t), since 
the growth with t is monotonic. 

We may factorise this result using eqn.(3.5): this gives 



J dr r"" 1 exp ( - r + dt'-y(N, a s (t + t'))) 



G(N,t). (3.12) 



The expression in square brackets is thus the resummed hard cross-section C(N;a s (t)). 
Note that if we Taylor expand G(N, t + r) in (3.11) in powers of r 



if 

n'.J a 



poo w m+n — 1 am 

dTT n - l e~ T G{N, t + r)= drV r-^e" T -z—G(N, t) 

y ' Jo n\m\ dt™ K ' 



m=0 

oo 



771 = 



(3.13) 



where in the second line we exchanged the order of summation and integration in order to 
perform the integrals over r. The result is of course precisely what one gets by first Taylor 
expanding 1/(1 — M) n in eqn.(3.9) about M = and then using eqn.(3.4). However the 
resulting series is useful only if 7(iV, t) and its derivatives are sufficiently small, as discussed 
in the previous section, whereas the integral representation eqn.(3.11) always converges. 
When the series diverges, the integral representation resums it. 

It is instructive to consider explicitly two analytic examples of how this works in 
practice. First consider what happens at fixed coupling: then at leading order G(N, t) = 
e y ^ N ' ,as ' t Go(N), so taking n = 1 for simplicity 

POO 1 

^ g (N,t)= / dTe- T e^ a '^G (N) = - -—G(N,t). (3.14) 

Jo l-7Wa s ) 

This is precisely the result expected from the pole dominance argument described in sec.2.4: 
M in the hard cross-section is replaced by 7(iV; a s ) according to eqn.(2.21). So the pole 
at M = 1 results in a new pole when 7(iV;a; s ) = 1, to the right of the rightmost pole in 
7(iV; a s ): for example if 7(iV; a s ) = a s /N, the new pole is at iV = a s . Thus on performing 
the iV integration (3.1) the cross-section E(p, Q) will exhibit a powerlike enhancement at 
high energy, growing faster than G{p, Q). 

Now consider what happens instead when the coupling runs: if we take for example 
1 {N;a s {t)) = a s (£)7o(iV), and a s (t) = l//3 t, then G(N,t) = t^ N ^^°G (N), and 

POO 

V* h (N,t) = / dTe-^t + ry o ^/^G (N) =t-^/ l3o e t T(l + ^^(^.^GiN.t). 
Jo 

(3.17) 
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io" 9 io" 8 io~ 7 io~ 6 io~ 5 io" 4 io~ 3 io~ 2 

Figure 9: a generic photo-production K- factor for a simple pole 1/(1 — M). The blue 
(upper) line corresponds to the exact expression eqn.(3.15), while the red, green and 
yellow lines (from bottom to top) are the NLO, NNLO and NNNLO approximation to 
it computed by including the second, third and fourth terms respectively of the series 
(3.16). The hard scale is set at lOGeV. 

where T(a, t) = J t °° dss a ~ 1 e~ s is the incomplete Gamma function. The only finite singu- 
larity in r(a, t) is a branch cut from t = down the negative real axis: for t real and 
positive, r(o,t) is entire in a. Consequently the only singularities of E;L(iV, t) are those 
of 7o(iV): in particular there is no new singularity when 7(£) = CK s (£)7o(iV) ~ 1. So when 
the coupling runs, S(p, Q) rises asymptotically in the same way as G(p, Q). 

It is easy to see that this smoothing away of the M = 1 pole (and indeed of any 
pole at M = mo, provided only that mo is real and positive) is a generic feature when 
the coupling runs: the key ingredient is that G(N,t) is a regular function of Int rather 
than of t. As recalled in sec. 2. 5 the interval M G (0, |) is then stretched to (0, oo), and 
so singularities at M = 1,2,... are effectively pushed out to infinity. This means that 
when the coupling runs the expansion of the hard cross-section around M = provides 
a good approximation to the exact result after only a few terms: the running coupling 
resummation of the M = \ singularity eqn.(2.29) in the evolution effectively deals with all 
the singularities for M > |, not only in the evolution but also in the hard cross-section, 
ensuring that all such singularities are factorised into the the initial distribution. 

In fact we can show that when the coupling runs, the expansion (3.13) is indeed 
asymptotic: returning to the specific example (3.17), and using the asymptotic series 

T(a,t) ~ t a-1 e~*[l + (a - l)t _1 + (a - l)(a - 2)t~ 2 + ...], (3.18) 
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as t — > oo, then as a 8 (t) — > 

^ h (iV, t) ~ [1 + 7 (t) + 7W 2 (1 - f ) + 7W 3 (1 - f )(1 - 2f ) + • • .]G(iV, f), (3.19) 

which is the same series term by term as (3.13) provided we evaluate the derivatives using 
eqns.(3.6). So in the running coupling case the expansion in powers of derivatives is indeed 
an asymptotic series. 

To see how well this asymptotic series works in practice, rather than in these simple 
examples, we define the If -factor as the ratio 1 

where £° h (p, Q) is simply the gluon distribution t). With this definition the de- 
pendence on the gluon distribution largely cancels, so one sees the effect of the hard 
cross-section. Then 

= 1 + cob s GK - 4) + cob S GK - t] + G<bS GK ' *> + ■ ■ ■ • <316 > 

where in the first line we used (3.11), and in the second (3.13). Using the resummed gluon 
distribution shown in fig. 4 and its derivatives fig. 8, it is now a simple matter to compute 
the if -factor, using either (3.15) or the series (3.16). The result is shown in fig. 9. It can 
be seen from the plot that even a simple pole at M = 1 gives rise to a quite substantial K- 
factor, and moreover that the NLO approximation to this if -factor, consisting of only the 
second term of (3.16), is adequate only for p ^ 10 -4 . However the convergence thereafter 
is very rapid: the NNLO approximation (the first three terms of (3.16)) is very good down 
to p ~ 10 -8 , and the NNNLO approximation (four terms) is difficult to distinguish from 
the exact result. Thus although the series is only asymptotic, in practice the first few 
terms give an excellent approximation to the full result. 

3. 3. Photoproduction of Heavy Quarks 

We now consider the inclusive cross-section for the photoproduction of a heavy quark 
pair. The hard scale in this case is the mass of the heavy quarks: Q = 2m q . The off-shell 
hard cross-section may be calculated from the diagram in fig. 5: the result in Q MS scheme 
is [35] 2 



nl AT 2 ,7V 14 + 20iV + 9N 2 + N 3 - M(10 + 7N + N 2 ) 

C(N, M) = enaa s 7r4 iV — — — — —7 

v ' Q 3-2M + 2N 

r(i- m + iv) 3 r(i + m) 



(3.21) 



x 



T(2 — 2M + 2iV)r(4 + N) 



Note that this if-factor is not the same as those used in phenomenological applications, 
which also include a factor due to the difference between LO and higher order partons. 
2 In refs. [15,16,35] this function is denoted by jwf/y). 



20 



When M = this reduces to the usual MS photoproduction cross-section relevant at large 
p, with poles at N = —1, — |, —2, . . ., while at N = it reduces to the impact factor [15] 

v ' ' Q 3 3-2M r(2-2M) v 7 

relevant for calculations at small p, i.e. at high energies S ^> Q 2 , with higher twist poles 
at M = - 1, — 2,... and infrared poles at M = 1, |, 2, . . .. It is the latter that are relevant 
for the high energy limit: in particular near M = 1 we have a double and a simple pole: 



2n 

C(0, M) ~ e|aa sy 7? -^ - + 0(1) . (3.23) 



2 1 

(i — m) 2 ~ r^M 



Expanding eqn.(3.22) about M = gives the Taylor expansion 

C(0, M) = 4aa.? (l + §M + ^ M 2 + - 2C 3 )M 3 + 0(M 4 )) 
~ e|aa s 2.444(1 + 1.952M + 3.87M 2 + 5.32M 3 + ...). 

This series has quite large coefficients, all of one sign, and growing: one thus expects it to 
converge rather slowly even for \M\ < 1. This is largely due to the double pole at M = 1: 
removing this by hand we find a series with rather smaller coefficients 

C(0, M) * e%aa s ^ (f ^ - § - §M - §M 2 + (lg - 2 C3 )M 3 + 0(M 4 )) 
~ egaa s 2.444( (1 1 : 7 ^ 4 )2 - 0.714 - 1.476M - 1.27M 2 - 1.54M 3 + ...). 

An even better series may be obtained by removing both the double and single poles by 
hand: 

C(0, M) ~ e%aa s 7 -f (fj^s - § ^ + \ - gM - §M 2 + (ff - 2C 3 )M 3 + 0(M 4 )) 

~ e|aa s 2.444( ^gp - jirg + 0.143 - 0.619M - 0.41M 2 - 0.68M 3 + ...). 

(3.26) 

We can now use these expressions to compute the if- factor for inclusive B photopro- 
duction, 

KM=b))' • (3.27) 



where S^ h (p, tub) is the cross-section computed using resummation, and S^ h (p, mg) the 
cross-section computed using the same gluon distribution, but with the LO hard cross- 
section (which here is simply C(0, 0) = ^-eQaa s , i.e. the first term in the expansion 
(3.24)). To calculate E^ h we use the same techniques as in the previous section: poles at 
M = 1 are dealt with using the exponentiation trick (3.11), while powers are turned into 
derivatives according to (3.4): the only new feature is that each term gets multiplied by the 
various coefficients in (3.24), (3.25), or (3.26). The results are plotted in fig. 10: the blue 
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Figure 10: The if-factor for photoproduction of b quarks. The blue curve (upper) is the 
resummed result, while the red, green and yellow curves (bottom to top) are the NLO, 
NNLO and NNNLO fixed order results computed as described in the text. The scale is 
Q = lOGeV. 



curves are the resummed results computed using the three approximations (3.24), (3.25), 
and (3.26). The results are indistinguishable on the plot, and thus may be taken to be the 
exact result: as in the previous calculation, the Taylor series is an adequate approximation 
to the more accurate pole approximations (3.25) and (3.26) for all p << 10 -10 . 

Also in fig. 10 we plot the if-factors for fixed order perturbation theory, at NLO, 
NNLO and NNNLO. Note that these are computed with the GLAP gluon distributions at 
the appropriate order (LO, NLO and NNLO), not with the resummed distribution. To do 
these calculations we use the result eqn.(2.27) to evaluate the dominant contributions to 
the fixed order hard cross-section: these are O(jf) at NLO (these are the terms computed 

in [36], not the full result of ref.[37]), O(^) and O(^f) at NNLO, and O(^), O(jfr), 

and O(^f) at NNNLO. They may thus be computed using the O(M), 0{M 2 ) and 0(M 3 ) 
terms in the Taylor expansion (3.24) provided we use at least the LO, NLO and NNLO 
results respectively for the evolved gluon distributions. It may be seen from the plot that 
while the NLO calculation seriously underestimates the X-factor for p ^ 10 -4 , the NNLO 
calculation is fine down to p ~ 10 -8 , while the NNNLO is worse, starting to overshoot 
below 10 -7 or so. So it seems that while a NLO calculation of the hard cross-section 
is inadequate at high energy, a NNLO calculation can perform well, the fully resummed 
result being only really necessary at very small values of p. 
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Figure 11: the resummation factor Rb for photoproduction. The solid blue (lower) 
curve includes resummation in both evolution and hard cross-section, while the dashed 
blue curve only has resummation in the cross-section. The green (upper) curve is the 
same computation comparing NNLO to NLO perturbation theory. 

To explore the relative sizes of the various resummation contributions, in fig. 11 we 
plot the "resummation factor" for photoproduction of bb pairs: 



S^u(p, TTIb) 
S^ u (p,m B ) 



where is computed as previously, while is computed using the NLO hard cross- 

section and NLO GLAP gluon distribution. The result is the solid blue curve: the net 
effect of the resummation is to reduce the cross-section by between 5 and 10%, rather 
uniformly over the full range of p. The dashed blue curve is the same calculation but 
with T,^ h computed using the NLO GLAP gluon, rather than the resummed one, to assess 
the relative effects of the resummation of the cross-section and the resummation of the 
evolved gluon: clearly both effects are of similar importance, since when combined they 
almost cancel. The green curve is the ratio £^ LO /£^°, with E^ LO computed with a 
consistent NNLO cross-section and gluon, for comparison with the resummed result: from 
this it is clear that a full NNLO calculation does not give a good approximation to the 
resummed cross-section for p ^ 10 -2 , despite giving a good account of the hard cross- 
section. This is because the NNLO evolution is not very close to the resummed evolution, 
as it underestimates the suppression [12]. Taken together, these three curves probably 
give a reasonable impression of the overall range of uncertainty in the 55-photoproduction 
cross-section at high energy. 
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4. Hadroproduction 

4-1. The Gluon-Gluon Luminosity 

When we calculate hadronic cross-sections, using the factorization (2.2) or equivalently 
(2.14), it is convenient to first define the gluon-gluon luminosity density 3 

L z (z, QuQi) = [ — / — S(z - x 1 x 2 )G(x 1 , Qi)G(x 2 , Q 2 ) 
J P xi J p x 2 

= J 1 ^-G(^,Q 1 )G(y,Q 2 ) (4.1) 

where £ = \nl/z, t x = lnQ 2 /A 2 , t 2 = \nQ 2 2 /p 2 . Thus we have in Mellin space simply 

L Z (N,M 1 ,M 2 ) = G(N,M 1 )G(N,M 2 ), (4.2) 

so eqn.(2.14) reads 

£ hh (p,Q)= r r d ^ d ^e^^H{N,M u M 2 )L z {N,M u M 2 ) 

J- ioo 2m J_ ioo 2m 2tvi 

(4.3) 

where now we have set t\ = t 2 = t = ln<5 2 /A 2 , Q being the invariant mass in the final 
state. Substituting (4.1) into (4.3) then gives back the gluon-gluon contribution to (2.2), 
as it should. 

It is straightforward to compute L z {p,Q,Q) = L z (p,Q) from the gluon distribution 
G(x 7 Q) shown in fig. 4: the result is show in fig. 12. Note that the luminosity only starts 
to rise at large Q when 10 or so, whereas for the gluon distribution the rise sets in 
already at x ~ 10 -3 . Again the resummed luminosity grows rather more slowly under the 
resummed evolution than with NLO GLAP evolution. 

Of course in a hadron collider it is not always possible to vary p and Q independently 
as one does in photoproduction: the inclusive cross-section for the hadroproduction of a 
final state of invariant mass Q depends instead on L Z (Q) = L Z (Q 2 /S, Q, Q), where S is the 
(fixed) centre-of-mass energy of the machine. This is plotted in fig. 13 for three different 
colliders: the Tevatron, the LHC, and a notional VLHC with ^f~S = 200TeV. Note that 
Q~ 1 L Z (Q) gives a rough estimate of the inclusive cross-section. 



In the literature it is more usual to call this quantity dC/dr, where r = X1X2: after integration 
C is then the total gluon-gluon luminosity. 
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Figure 12: The gluon-gluon luminosity density L z (p,Q) plotted against Q in GeV for 
p = 1CF 10 , 1CP 9 , . . . , 0.01, 0.1 (from top to bottom). The blue curves are evolved with 

the NLO resummation described in the text: the black and red curves are with LO and 
NLO GLAP evolution. 
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Q/GeV 

Figure 13: The gluon-gluon luminosity L Z (Q) plotted against the invariant mass Q in 
GeV. The lower curves are for the Tevatron, the middle curves for LHC, and the upper 
for a VLHC: the colour coding is the same as in fig. 12. 

4-2. Resumming the Dominant Infrared Singularity 

In order to compute the inclusive cross-section eq.(4.3) we must first compute the 
hard cross-section, and then convolute it with the gluon-gluon luminosity. In practice 
this means we must perform the inverse Mellin transforms over N, M\ and M 2 . Since all 
the collinear and high energy logarithms are already included in the luminosity, the hard 
cross-section H(N, Mi, M2) is regular at N = Mi = M2 = 0, and we may Taylor expand 
it using (2.24). Just as in the photoproduction case eqn(3.4), the resulting integrals over 
powers of Mi and M2 may then be evaluated in terms of derivatives of the luminosity: 
writing L Z (N, t) = L Z (N, t, t) for simplicity, then for all mi, mi = 0, 1, 2, . . . 



dM 1 dM 2 
2-ni 2m 



e t(M 1+ M 2 ) M pMpL z (N, Mi, M 2 ) 



Qm 1 +m 2 

dt^dt™ 2 



L z (N,ti,t 2 ] 



ti=t 2 =* 

(4.4) 
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Figure 14: The first (left) and second (right) derivatives with respect to t = lnQ 2 /A 2 of 
the gluon-gluon luminosity L z (p, Q) plotted in fig. 13, normalised to L z , plotted against 

Q in GeV, for p = 10~ 10 , 1CF 9 , . . . 0.01 (from top to bottom). Again the blue curves are 
with NLO resummed evolution, the red with NLO GLAP evolution. 



The derivatives of the luminosity are then given in turn by 



d_ 

dfx 
d_ 2 

dt'f 



L Z {NMM) 



-^L z (N,t u t 2 ) 



U=u=t 



_d_ 

df 2 
d 2 



L Z (N, ti, t 2 ) 



-^L z (N,t 1: t 2 ) 



a 2 



a 3 
dtf 



-^L z (N,h,t 2 ) 



dt x dt 2 
d 3 



L g {N,t u t 2 ) 



t 1= t 2 =t dto 



d 3 

^rL z {NMM) 



dt(dt 2 ' 



-7^L z {N,t u t 2 ) 
t!=t 2 =t dt\dt\ z ^ ' 1 ^ 2 ' 



h=u=t 



u=t 2 =t 



u=u=t 



h=u=t 



u=u=t 



l{N,t)L z {N,t), 
(j 2 + j)L z (N,t), 
^L z (N,t), 

(7 3 + 3 7 7 + 7)^(^,t), 
( 7 3 + 77 )L 2 (iV,t),... 



(4.5) 

The first two of these are plotted in fig. 14 (to be compared with fig. 8): again even the 
first derivative is below unity for all except the smallest values of p and Q, and the second 
derivative is considerably smaller than the first. The resummed derivatives are smaller at 
small p and Q than the corresponding GLAP derivatives. Given these facts, we expect the 
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Figure 15: A schematic picture of the location of the singularities of a hadroproduction 
impact factor H(0, Mi, M2) in the M1-M2 plane. The particular singularities shown are 
those for heavy quark production: those in green (parallel to the axes) are simple poles, 
in blue double poles and in red triple poles (diagonal). 

first few terms of the Taylor series to give a good approximation to the full cross-section 
at all except the highest energies and lowest scales, just as they did in photoproduction. 

It remains to consider the singularities of the hard cross-section (see fig. 15).. Just 
as in the photoproduction case there will be higher twist singularities at Mi,M2 = 
— 1, —2, —3, . . . which will lead to terms suppressed by inverse powers of the hard scale 
Q 2 : these will not concern us here. The structure of the infrared singularities is more 
subtle: there are singularities at Mi, M 2 = 1, 2, 3, . . . just as in photoproduction, but there 
are now also infrared singularities on the lines Mi +M2 = 1, 2, 3 . . ., as discussed in sec. 2. 6. 
Since these come closest to the origin, it is likely that they dominate at high energies, and 
indeed as we shall see this turns out to be the case. 

In order to integrate over the line of poles at Mi + M2 = 1, we may employ the same 
trick that we used for the M = 1 pole in sec. 3. 2: 



for n = 1, 2, The nice extra feature here is that under the integral over r the dependence 




(4.6) 
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on Mi and M 2 has now factorised, allowing both integrals to be performed independently: 



ZUN,t) = f°° ^Ml^Ml 1 -e t ^ +M ^G(N,M 1 )G(N,M 2 ) 

hhV ' ; J_ ioc 2m 2m (1 - M x - M 2 ) n V ' V ' 

= J_ r dTT n-l e -r f l0 ° ^Ml e (t + T)M lG{ ^ Mi) ^(t+^M^^^ 

1 Z" 00 

= / drr^e-^L^t + r), 
ni Jo 

(4-7) 

where in the last line we have exploited the factorization to perform both of the integrals 
over Mi and M 2 , writing the result in terms of the gluon-gluon luminosity L Z (N, t). In this 
way we reduce a double integral over two complex variables with a line of singularity to 
a single real integral, which clearly converges rather rapidly when the coupling runs since 
L z (N,t) is very smooth (see fig. 12). Again where L z (p,Q) increases monotonically in Q 
the singularity will give an enhancement of the cross-section, which will be most significant 
when the rise is steepest (i.e. small p and low Q). 

It is easy to see that the same arguments used for the photopro duct ion and electro- 
production cross-sections, in particular (3.14) and (3.17), will also apply here, because 
of the factorization under the r integral in (4.7). Thus when the coupling runs the line 
singularity is once again pushed to infinity, and the only singularities of £^ h (iV, t) in the 
iV-plane are those of the anomalous dimension 7(iV, t). The running of the coupling thus 
ensures that all infrared singularities, whether in the evolution or in the hard cross-section, 
are factorised into the initial (nonperturbative) gluon distribution. 

An important consequence of this factorization is that the asymptotic high energy 
behaviour of inclusive hadroproduction cross-sections is determined entirely by the rise in 
the gluon distributions, just as it was in electroproduction and photoproduction. It follows 
that the high energy powerlike rise due to the rightmost singularity in the anomalous 
dimension is universal: all inclusive cross-sections rise in the same way. 

To see how well all this works in practice, consider a triple pole, i.e. n = 3, and define 
the X-factor as the ratio 

where the denominator E° /l (Q 2 /5', Q) is simply the gluon-gluon luminosity distribution 
L z (Z,t) when f = InS/Q 2 . Then 

1 r io ° rlN f°° 

K ^Q) = -TT^ 7T- eNC drr 2 e- T L z (N,t + r) (4.9) 



L z (£,t) J- iQ0 2ni J 

3 d T _ , 6 d 2 T _ , 10 d 3 T _ , 

- 1 + mj)Wt L ^ + TMJ)W L ^ + mJW^ + " <410) 

where in the first line we used (4.7), and in the second we Taylor expand L Z (N, t + r) and 
integrate term by term. The diagonal derivatives of the luminosity are readily deduced 
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Figure 16: a generic hadropro duct ion X-factor for a triple pole 1/(1 — Mi —M2) 3 , plotted 
against invariant mass Q, for three different colliders: the Tevatron (lower curves), the 
LHC (middle curves) and a VLHC (upper curves). In each case the blue (upper) lines 
correspond to the exact expression eqn.(4.9), while the red, green and yellow (bottom to 
top) lines are the NLO, NNLO and NNNLO approximation to it computed by including 
the second, third and fourth terms respectively of the series (4.10). 



from (4.5): 



a 2 

dt 



2 L z (N,t) 



Q3 



L z [N,t) 



2-y{N,t)L x {N,t), 
2(2 7 2 + 7 )L 2 (iV,t), 
2(4 7 3 + 6 77 + 7 )L z (iV,t), 



(4.11) 



Using the resummed luminosity distribution shown in fig. 12 it is now a simple matter to 
compute the .fT-factor, using either (4.9) or the series (4.10). The results are shown in 
fig. 16, for S = 2 TeV, 14 TeV and 200 TeV. At sufficiently low invariant mass Q and thus 
large p the K factor rises steeply as expected, eventually resulting in very large corrections. 
In each case the exact result (4.9) is well approximated by only the first few terms in (4.10). 
In fact above Q ~ 5 GeV at the Tevatron, Q ~ 15 GeV at the LHC and Q ~ 70 GeV at 
a VLHC the O(M) (ie NLO) correction alone already gives a good approximation: only 
below these scales do the 0(M 2 ) (NNLO) corrections become significant. The NNLO 
approximation only really starts to be inadequate when Q 5 GeV at the VLHC, which 
would no doubt be outside the range of acceptance. 
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Figure 17: the same as fig. 16, except that now the red, green and yellow lines are the 
NLO, NNLO and NNNLO results in fixed order perturbation theory (at high energy), 
computed by including the second, third and fourth terms respectively of the series 
(4.10), with the luminosity evolved using GLAP evolution at the appropriate order. 



These results are interesting because, as shall show explicitly for heavy quark produc- 
tion in the next section, the line of triple poles (2.32) is sufficiently singular to provide the 
dominant contribution to the cross-section. Moreover since these infrared singularities are 
quite generic, we expect other inclusive cross-sections such as the inclusive jet cross-section 
to acquire iiT-factors at high energy similar to those shown in fig. 16. Calculation of the 
.KT-factor can then proceed in each case through a computation of the expansion of the hard 
cross-section H(N, Mi, M2) in Mi and M2, ie by perturbation about the on-shell result: 
the singularity will then generate large coefficients of the 0(Mi + M 2 ) and 0((Mi +M 2 ) 2 ) 
terms, and thus large NLO and (at high enough energy) NNLO X-factors. 

This mechanism thus provides a simple explanation for the large X-factors commonly 
found in hadroproduction processes at high energy. Constraining the incoming gluons 
to be on-shell as in the usual LO calculation keeps the timelike intermediate gluon in 
fig. 6b or fig. 7b away from its mass-shell. Releasing this constraint either by using the 
off-shell formalism used here, or by going to higher order in a s in the more usual on-shell 
formalism (with in particular contributions from diagrams in which one (NLO) or both 
(NNLO) incoming gluons emits another gluon), allows the intermediate gluon to get close 
to its mass-shell, and thus produces the large enhancements evident in fig. 16. 

To test this idea it is instructive to compare the resummed X-factors with those 
computed in fixed order-perturbation theory. At high energy this is easily done by using 



31 




Figure 18: the resummation factor R3 for a generic hadroproduction process. As in 
fig. 11 the solid blue (lower) curve includes resummation in both evolution and hard 
cross-section, while the dashed blue curve only has resummation in the cross-section. 
The green (upper) curve compares NNLO to NLO unresummed perturbation theory. 



the expansion (4.10), applied to the luminosity evolved using GLAP evolution (at LO for 
the NLO if-factor, NLO for the NNLO, and NNLO for the NNNLO), just as we did for 
photoproduction in fig. 10: the results are shown in fig. 17. At the Tevatron NLO GLAP 
perturbation theory for the hard cross-section does pretty well except at very low scales 
(appropriate for charm production). At LHC NNLO is starting to be important already 
around the beauty threshold, while at a VLHC NLO is only good at scales above 60 GeV. 
For beauty production at VLHC the NNLO correction is very large - around 50% above 
NLO. More importantly, the NNNLO correction is nowhere an improvement on the NNLO 
result: beyond NNLO the resummed result is more useful. 

To assess the impact of the resummation, in fig. 18 we plot the "resummation factor" 

T>h(ri\ — hh 

where, as in eqn.(3.28), the reference cross-section Sj^^* is computed using the NLO hard 
cross-section and NLO GLAP gluon distribution. Again we have a substantial cancellation 
between the suppression of the luminosity due to the resummation of the evolution and 
the enhancement of the cross-section due to the triple pole in the hard cross-section. At 
intermediate scales (Q ~ 8, 20, 60 GeV for Tevatron, LHC and VLHC respectively) the 
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resummation gives an overall suppression of around 10%: below these scales the suppres- 
sion starts to turn into an enhancement due to the rapidly rising partonic cross-section. 
However these enhancements are generally at scales sufficiently small that they would es- 
cape detection in all but a very forward detector. Again we also show the ratio of the 
NNLO to NLO fixed order computation: this overestimates the cross-section because it 
underestimates the suppression of the gluon-gluon luminosity due to resummation. 

4-3. Hadroproduction of Heavy Quarks 

As a specific example of resummation in hadroproduction we now consider the 
hadroproduction of heavy quarks, specifically 66-production at the Tevatron, LHC and 
VLHC. The LO contribution to the hard cross-section then involves the computation of 
the two diagrams in fig. 6 with both incoming gluons taken off-shell. Unlike the photo- 
production cross-section eqn.(3.21) the hadroproduction cross-section has been evaluated 
analytically only in the high energy limit A = 0: the result is the impact factor [16,27] 4 

h(o,m u m 2 ) = a 2_^_r(i + Ml) r(i _ m 1 )t(i + M 2 )r(i - m 2 ) 



X 



f (r(3- M 1 -M 2 )) 2 / / r(l-Mi-M 2 ) \ 

c (l - Mi - M 2 )r(6 - 2(Mi + M 2 )) v \r(i -Mi)r(l -M 2 )J 



A c v y ' ' r(4 - 2Mi)r(4 - 2M 2 ) 

The second term in this expression is due solely to the "abelian" diagram fig. 6a: when one 
leg is on-shell (e.g. M 2 = and M\ = M) this term reduces to C(0, M) eqn.(3.22) up to an 
overall constant vertex factor. The first term is due to the intrinsically nonabelian diagram 
fig. 6b: it is this piece which contains the triple pole singularity (2.32), which dominates 
the cross-section at high energy. 

Consider the structure of the result eqn.(4.13) in various regions of the M\-M 2 plane 
fig. 15. It has higher twist (simple) poles at Mi,M 2 = —1,-2,..., and infrared (anti- 
collinear) poles at Mi, M 2 = 1,2,...: for example near Mi = 1 

HftM,,^)-^! 1 r(M a )r(i-M,)(r(2-M a) )' i 

V ' U 2) s A c Ll-M! T(4-2M 2 ) v } y v ; 

thus a simple pole except when M 2 is an integer, at which special points there is a double 
pole. However it also has lines of singularity when Mi + M 2 = 1,2, . . .: writing M± = 
Mi ± M 2 , the Laurent expansion about M + = 1 is 



H(0,M u M 2 )^a 2 s j^- I 



N c I -Ml 



6 (1-M+) 3 

+ jr^JurTi( 3{1 - M -^v - h\ + \m_) - m - \ M _)) - 11 + 5Mi) 

- ( (67+72(lD2) ;- 132lD2)JVc - 3^: + 0{M1)) + 0(1) 



(4.15) 



In refs. [15,16,27] the function H(N,M 1 ,M 2 ) is denoted by ^(71,72). 
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i.e. triple, double and simple poles, except again at the special points M_ = ±1 (and thus 
(Mi, M 2 ) = (1, 0) or (0, 1)) where the triple pole reduces to a double pole, the double to a 
single. 

For Mi and M 2 close to zero the impact factor is regular, as it must be, with Taylor 
expansion 

ff(0, Mi, M 2 ) = al-^- [(^ - * ) + - « ) (Mi + M 2 ) 



, / 4924AT e _ ^22_Vm 2 + M 2 1 + ( (9848-150tt 2 ) 7V C _ (470+217T 2 ) \ , f , f 
+ 1 3375 81N C )\ m l + 1V1 2> + I 3375 324N C ) m l m 2 



2 



(4.16) 



, ^15U544-27UUU<;3jiVc (730-189(,3) \ / A/f 3 , A/f 
+ I 50625 2431V; ) + M : 

+ ( (150544-22500C3-1925^)iV c _ (2776+378^+123^) ) (Mj + M2)MlM2 + 0(M 4 } 

~0.2633a 2 (l + 2.69(Mi + M 2 ) + 5.78(Mi + M 2 ) 2 - 1.50MiM 2 

+ 9.41(Mi + M 2 ) 3 - 2.79(Mi + M 2 )M 1 M 2 + 0(M 4 )). 

The rather large numerical coefficients in this expansion (compare them to those of the 
expansion (3.24)), in particular of the powers of Mi +M 2 , are due to the dominance of the 
nearby triple pole singularity (4.15): if we first subtract the triple pole, and then Taylor 
expand what is left, we find 

H( ,M u M,),a* 7 ^- 1 [^ {1 _ M 1 i _ M2r + (ft - 

+ " + M 2 ) + (if^ - ^)(M? + M|) 

, / (3098-1507r 2 )jV c _ (470+21tt 2 ) \ M M 
I 3375 324Af c ) lvl l lvl 2 

, / (66169-27000(3)^ (730-189^) W A/f 3 , , f 3\ (4.17) 
+ I 50625 243A^ j + iw 2 J v ' 

+ ( (66169-22500C3-1925^)AT C _ (2776+37^+123^) ) ^ + + 0{M*)] 

-0.2633a 2 ( (1 _° A J i 4 ^ 2) 3 + 0.2541 + 0.448(Mi + M 2 ) + 1.31(Mi + M 2 ) 2 

- 1.50MiM 2 + 1.95(Mi + M 2 ) 3 - 2.79(Mi + M 2 )MiM 2 + (9(M 4 )) . 

A further reduction in the coefficients may be obtained by also subtracting the double 
pole. 

It is now straightforward to compute the hadronic cross-sections using either the 
Taylor expansion (4.16) or the more precise expansion (4.17), the expressions (4.4) and 
(4.5) for the powers of Mi and M 2 , the representation (4.7) for the triple pole in (4.17), 
and the gluon-gluon luminosity shown in fig. 12. The results are presented in fig. 19 as a 
K- factor 

KBM =mgm„ .... (^8) 



Q 2 =m] l+p 2 



where £jf h is the fully resummed calculation, and the reference cross-section TP hh is eval- 
uated using i?(0, 0,0) = a^MM' tnus dividing out the overall normalization and the pri- 
mary dependence on the gluon-gluon luminosity. Both cross-sections are evaluated at 
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Figure 19: the i\-factor for hadroproduction of bb pairs as a function of their minimum 
Pt, at the Tevatron (lower), LHC (middle) and VLHC (upper). As in fig. 17 the blue 
curves are the resummed result, while the red green and yellow curves are fixed order 
perturbation theory at NLO, NNLO and NNNLO respectively, computed as described 
in the text. 

Q 2 = m 2 B + p T , where px is its minimum average transverse momentum of the bb pair. 
Note that we do not compute the transverse momentum distribution (the hard cross-section 
H(N, Mi, M%) is fully inclusive): the pr dependence of (4.18) is simply a reflection of the 
restriction of the phase space when we require a higher invariant mass in the final state. 
We expect this to be the dominant effect here. 

As we found in the previous section the naive Taylor expansion (4.16) and the pole 
resummed expansion (4.17) give almost identical results throughout the entire kinematic 
range: the two resummed curves in fig. 19 are indistinguishable. In fact the convergence 
is so rapid that only the 0(M) (NLO) and 0(M 2 ) (NNLO) terms in (4.16) are actually 
needed: NNNLO contributions would only become significant for charm at a VLHC. 

For comparison we also show in fig. 19 results for fixed order (GLAP) perturbation 
theory at NLO, NNLO, and NNNLO, computed using the high energy approximation 
(as in fig. 17) by using GLAP evolution, and keeping only the O(M), 0(M 2 ) and 0(M 3 ) 
terms respectively in (4.16). At the Tevatron both the resummation and the NNLO GLAP 
corrections to the partonic cross-section give a slight suppression, while at the LHC there 
is a modest enhancement, at least at low px- Only at the VLHC is there a substantial 
enhancement, and even there a fixed order NNLO calculation of the partonic cross-section 
would be quite sufficient. 

In order to estimate the overall effect of the resummation, we show in fig. 20 the 
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Figure 20: the resummation factor Rb for a hadroproduction of bb pairs plotted against 
their transverse momentum. As in fig. 18 the solid blue (lower) curve includes resum- 
mation in both evolution and hard cross-section, while the dashed blue curve only has 
resummation in the cross-section. The green (upper) curve is the same computation 
comparing NNLO to NLO perturbation theory. 



resummation factor 

^ h (Q 2 /S,Q) 



T$Z°(CP/S,Q) 



Rb(Pt) = u'a^ ^ 2 ^ 2 - ( 4 - 19 ) 



Q 2 =m 2 +p 



The reference cross-section is now the NLO fixed order cross-section computed with the 
NLO GLAP evolved gluon-gluon luminosity: unlike in the X-factor effects due to the 
different evolution of the luminosity are thus now included. As in fig. 18 we find that the 
overall effect of the resummation is a net suppression by about 10%, which gradually goes 
away as pr increases. Just above threshold this suppression is remarkably independent of 
the centre-of-mass energy S of the machine, and when this is very high (in particular for 
VLHC) the suppression goes away very slowly. So once again the enhancement expected 
from the triple pole in the cross-section is more than compensated by a suppression of the 
gluon-gluon luminosity. 

It should be noted that all these calculations are only estimates: in particular a 
proper matching to the high p (Sudakov) region, inclusion of quark effects and realistic 
fitted parton distributions could all change the results substantially. In particular the 
suppression of the gluon-gluon luminosity by the resummation of the evolution is probably 
overestimated, since the starting distribution at 2 GeV is held fixed, rather than fitted 
to data. However the band between the upper and lower curves in fig. 18 and fig. 20 is 
probably a reasonable estimate of the current overall uncertainty due to resummation in 
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hadroproduction cross-sections: for inclusive B production this means roughly —5 ± 5% 
at the Tevatron, 5 ± 10% at the LHC and 20 ± 20% at a VLHC. In all cases it is probably 
comparable to the impact of NNLO in the partonic cross-section, and may result in a net 
suppression rather than an enhancement. At large pr the effect goes away as expected. 

4-4- Drell-Yan and Higher Orders 

Of course not all processes have the same structure of infrared singularities as heavy 
quark and inclusive jet production: as the number of particles in the final state increases, 
so does the range of possible infrared singularities. Indeed, adding one more particle to the 
final state general adds one extra collinear and one extra soft singularity, thus increasing 
n in eqns.(4.6) and (4.7) by two. 

Consider for example the gluonic contribution to Drell-Yan or vector boson produc- 
tion, given by the diagrams in fig. 3. The most singular of these diagrams is fig3c: there is a 
soft singularity from the timelike gluon, a collinear singularity from the splitting into a qq 
pair, and then further soft and collinear singularities from the final vector boson emission. 
Unless there are accidental cancellations, we thus expect the impact factor for this process 
to have a pole of order five at Mi +M 2 = 1. The relevant i^-factor is thus (using eqn.(4.7) 
with n = 5) 



which may be evaluated either by computing the r integral numerically, just as we did for 
K 3 eqn.(4.9). 

The results are shown in fig.21. Comparing K§ to K% shown in fig. 17, the extra soft 
and collinear singularities produce a further overall enhancement as expected. However 
the qualitative features of the two plots are very similar: in particular the resummed result 
starts to grow faster than the NLO result at about the same scale in each case. For W 
or Z production at LHC the correction is in the region of 20-30%: it only becomes large 
at VLHC. However for production of Drell-Yan pairs at around lOGeV at LHC the NLO 
correction is as large as a factor of three, and requires resummation. 

For these gluon-gluon processes NLO means O(aal), which is NNNLO in the usual 
nomenclature of fixed order perturbation theory. Thus unlike in the previous heavy quark 
and inclusive jet (ie n = 3) estimates, here even the NLO curve is a new result: only the 
LO contribution to the gluonic contribution to Drell-Yan and vector boson production, 
gg — > W + X (i.e. the graphs in fig. 3 but with the incoming gluons on-shell), has been 
computed exactly in fixed order perturbation theory [20]. 

We can also use fig.21 to estimate resummation corrections to prompt photon pro- 
duction, or to the three-jet inclusive cross-section, since the relevant diagrams again have 
the same structure as those in fig. 3 (with for jets the quarks and vector bosons replaced 
by gluons) and thus the same structure of infrared singularities. 

Since the number of infrared logarithms increases by two at each extra order in a 8 , and 
since, unlike the logarithms from high energy (N = 0) and collinear (M = 0) singularities 
these logarithms are not being explicitly resummed, one might worry that they might 




(4.20) 
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Figure 21: the Zf-factor K5, eqn.(4.20), appropriate for the gluonic contribution to Drell- 
Yan, vector boson production and prompt photon processes. The curves are labelled as 
in fig. 17: blue is resummed, while red and green are estimates of fixed order perturbation 
theory at NLO and NNLO. 



destabilise the hierarchy of terms in the resummed perturbation theory [15]. To address 
this problem we note that we can estimate the size of subleading resummation corrections to 
by computing the ratios a s (Q)K n+ 2(Q) / K n (Q): these are plotted in fig. 22 for n = 3 (heavy 
quark production and inclusive jets) and n = 5 (Drell-Yan and vector boson production). 
It is clear from this plot that even at low scales and at VLHC the enhancements due to 
the two extra infrared logarithms are never very substantial. 

Moreover as n increases the enhancements are systematically reduced. The reason 
for this is not hard to find: a singularity of order n results in a smearing of the gluon- 
gluon luminosity with a distribution proportional to r n_1 e _T , eqn.(4.7), which is peaked 
at r = n — 1. Since the scale dependence of the luminosity is very smooth, thanks to 
asymptotic freedom, the main effect of the smearing at large n is to shift the scale at 
which the luminosity is evaluated from t to t + n — 1 . Thus at large n 

K„+2(Q) LM,t + n+\) 2 B 

-kJW £,«, t + » - 1) ~ 1 + m, < + ») a i £,({ ' ( ' 

which tends to one for large n and large t (see fig. 14). It follows that there is no reason 
to suspect that the infrared logarithms spoil our resummed perturbation theory: for high 
energy processes with a single hard scale, we have indeed resummed all large logarithms. 
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10 100 1000 

Q/GeV 

Figure 22: the ratios a s K n+ 2/ K n for n = 3 (blue, upper) and n = 5 (green, lower) 
showing relative enhancement of subleading corrections. As usual the three curves of 
each colour are for Tevatron, LHC, and VLHC. Also shown is a s {Q) (solid black) as a 
baseline expectation. 

5. Rapidity Distributions 

5.1. Gluon-gluon Rapidity 

Besides hadroproduction total cross-sections it is also interesting to consider rapidity 
distributions: in the central rapidity region both partons carry roughly the same fraction 
of longitudinal momentum, but at large rapidities one of the partons is at a much smaller 
value of x than the other (see fig. 1 for the ranges covered at various machines) so it is 
perhaps here that one might expect the effects of resummation to be most significant. 

We define 

z = x\X2, T] = \ ln(xi/x 2 ), (5.1) 

so that s = zS is the centre-of-mass energy in the partonic collision, and rj is the (pseudo)- 
rapidity: 77 = in the central region, becoming large and positive/negative in the for- 
ward/backward regions. In terms of z and 77 the fraction of longitudinal momentum in 
each of the two colliding partons is 

xi = v^e 77 , x 2 = v^e" 77 . (5.2) 

Since s > Q 2 , while xi, x% < 1, we must have 

p<z<l, ilnz < ry < |lni, (5.3) 
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Figure 23: The gluon-gluon rapidity distribution L^(ry, Q) at the LHC for 5-production 
(Q = lOGeV) (upper curves) and VF-production (Q = 75GeV) (lower curves). The blue 
curves are for gluons evolved using NLO resummation, while the black, red and green 
curves are for gluons evolved using LO,NLO and NNLO GLAP. 

In terms of z and r\ the factorization formula (2.2) may be written 

Shh( p,Q) = f p d ij}% /^/^(e, 1^)0(^,^(^,1^). 

(5.4) 

Since the hard cross-section depends only on z, not on rj (because of invariance under 
longitudinal boosts), we may perform the integral over rj first to give the gluon-gluon 
luminosity (4.1). If instead however we perform the integrals in the opposite order 

" \ drj= f drj f — , (5.5) 



z 



2 inz 



where as usual £ = lnl/p, the factorization (5.4) may be written instead as a factorization 
for the differential cross-section: 



az^hh 



drj 



d i 7$/S s » ( ^ > l )G( ^' k?)G( ^" , ' ki) - (5 - 6) 
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At high energy the explicit dependence of the hard cross-section £ 3S (f , ^) on z 
is relatively weak: the dominant contribution to the z dependence of the collinear cross- 
section comes about through the dependence of the off-shell cross-section on the transverse 
momenta. This makes the computation of rapidity distributions at high energy particularly 
simple, since the rapidities of the final state particles are directly related to the rapidities 
of the colliding partons. It is then useful to define gluon-gluon rapidity distribution 

L v ( V ,k 2 ,k 2 ,p) = / -G{V~ze\kl)G{V~ze-\kl). (5.7) 
J p z 

Clearly the normalization of this distribution is not independent of that of the gluon-gluon 
luminosity (4.1): in fact 

/ d V L v ( V ,klklp)= — L z (z,k 2 ,k 2 ). (5.8) 
J-ZI1 J P z 

Again the gluon-gluon rapidity distribution may be readily computed from G(x, Q) dis- 
played in fig. 4: the result for L v (r],Q) = L v (r], Q 2 , Q 2 , p) is shown in fig. 23 for both 
£?-production and VF-production at LHC. The falloff as r\ — > ±£/2 is very rapid, essen- 
tially because the gluon is very small at large x (remember that in all these calculation we 
are suppressing the quark contribution, so the valence region is underpopulated). 

In terms of the gluon-gluon rapidity at high energy the differential cross-section is 
then given by 



cEhh fd 2 \s.i fd 2 k 2 ^ ki k 2 



^kf J ^ E "^~Q'~Q^ L ^ k ^ k2 ' p) 



(5.9) 



2m 2ni 

where the (double) Mellin transform is defined in the usual way: 

f°° elk 2 /k 2 \~M 1 r°°fJk 2 /k 2 \~M 2 

L v (r,,Mi,M 2 , P ) = jf d ^-{%) I ^(g) L V ( V MM,P). (5.10) 

Note that unlike (4.3) the expression (5.9) only holds in the approximation where we ignore 
the N dependence in the hard cross-section H(N, M 1; M 2 ). 



5.2. B -production and W -production at LHC 

For a particular impact factor H(0, Mi, M 2 ) the corresponding differential cross- 
section (5.9) may be readily evaluated using eqns.(4.5) for powers of Mi and M 2 , and the 
exponentiation trick eqn.(4.6) for singularities. We consider two examples: S-production 
at the LHC, for which we use the impact factor (4.13) expanded as (4.16) or (4.17), and 
the gluon-gluon contribution to VF-production at LHC, which we assume is dominated at 
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Figure 24: the if -factors K Tt { r q,mB) (upper curves) K v (r),mw) (lower curves) for the 

rapidity distribution of -B-production and (the gluon-gluon component of) VF-production 

at the LHC. As in fig. 17 and fig. 19 the blue curves are the resummed result, while the 

red and green curves are NLO and NNLO perturbation theory, estimated in the usual 
way. 



high energy by the infrared singularities in final state emission (fig. 3c), and may thus be 
modelled by an n = 5 pole: 



jyW 

dq 



aa 2 a r 



lOO 

oo 



dMi dM 2 c t (Ml +M 2 ) 
2ni 2ni 



1-Mo-Mo) 5 



L v ( v ,M 1 ,M 2 ,p) 



Q=m\ 



p oo 

/ dTT A e- T L v (r),t + T) 
Jo 



(5.11) 



t=\nm 2 w 



/A 2 



where r is a normalization factor (just a number), and in the second line we have used 
eqn.(4.6) to evaluate the integrals over Mi and M 2 just as we did in eqn.(4.20). For B- 
production we use a similar expression, but here of course we take n = 3. Note however 
that now the integration over z in eqn.(5.7) can take us into the region of large X\ or x 2 
at high rapidities, so here we must take care to match the cross-section smoothly in this 
region to the fixed order calculation to avoid spurious contributions. 

Again we express the results of these calculations as if-factors: here 



K v (r] : m w ] 



a2j hh 



drj I drj 



Q=m w 



(5.12) 
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Figure 25: the resummation factor R v for ^-production (outer curves) and W- 
production (inner curves) at the LHC plotted against the rapidity rj. The solid blue 
(lower) curves correspond to fully resummed calculations, while the dashed curves are 
the same calculations but with the effect of the evolution on the gluon-gluon rapidity 
distribution factored out. The true size of the resummation effect probably lies some- 
where between these two extremes. The green (upper) curves are the result of a NNLO 
fixed order calculation. 

and a similar expression for K^(rj, mB)- The denominator is in each case the result ob- 
tained by setting Mi = M2 = in the impact factor, and thus divides out the unknown 
normalization r and the primary dependence on the gluon-gluon rapidity. The correspond- 
ing iiT-factors in perturbation theory may be estimated by expansion in powers of Mi and 
M2, and using the appropriate GLAP evolved gluon-gluon rapidity distribution. 

The results of these calculations are shown in fig. 24. For S-production the X-factor 
is already substantial in the central rapidity region, and rises further at large rapidities 
where one of X\ or £2 becomes very small. This rise is rather steeper for the resummed 
and NNLO calculations than it is at NLO, as expected. However this prediction needs 
to be interpreted with care, since at large rapidities one of the gluons is moving towards 
the high x region, which is poorly modelled in this calculation since there are no quark 
contributions. For the gluon-gluon contribution to W production the i^-factors in the 
central region are all rather small, but increase quite quickly with rapidity. However here 
the resummed calculation is close to the NLO and NNLO calculations throughout the 
whole range. Note again that here even the NLO curve is a new result. The change in 
the shape of the rapidity distribution due to what is formally a NNNLO (i.e. 0(ao? 8 )) 
contribution is thus quite striking. 
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To see more clearly the effect of the resummation alone, we also compute the "resum- 
mation factor" 



again with a similarly expression for 5-production, where now the reference differential 
cross-section is the NLO fixed order cross-section computed with a NLO GLAP evolved 
rapidity distribution, so that the effect of the resummation in both evolution and partonic 
cross-section are combined. The results are shown in fig. 25. Resummation reduces the B 
cross-section by around 10% across the whole rapidity region, this effect being due almost 
entirely to evolution since without it there is an enhancement rising to around 15% at large 
rapidity. The effect of resummation on VF-production is rather less pronounced. Note that 
in both cases the enhancement due to the infrared singularities in the hard cross-section is 
largely cancelled by the suppression of the luminosity, leaving a relatively flat distribution. 



6. Conclusions 

We have shown that the problem of integration of infrared singularities in photo- 
production, electroproduction and hadroproduction cross-sections due to final state gluons 
becoming soft or collinear [15,16] may be solved through an exponentiation trick eqn.(3.11) 
and eqn.(4.6) respectively. This enables us to show that despite the dramatic enhance- 
ments found when when the coupling is fixed, when the coupling runs the effects are much 
more modest. In particular we have shown that when the coupling runs the growth of the 
inclusive cross-sections at asymptotically high energy is given universally by the growth of 
the resummed integrated gluon distribution, with no process dependent power enhance- 
ments: all the infrared singularities are correctly factorised into the initial nonperturbative 
distribution, as they should be. 

We find furthermore that except at very low scales the inclusive cross-sections are well 
approximated by keeping only the first few terms in the Taylor expansion of the partonic 
cross-section in powers of M, irrespective of the presence of the nearby infrared singu- 
larities. Since this expansion in powers of M is closely related to the usual perturbative 
expansion of the hard cross-section in powers of a s /N, this enables us to understand the 
behaviour of the expansion to fixed orders in a s . In particular we have shown that the large 
.KT-factors commonly found in hadroproduction processes at NLO and NNLO are at high 
energy due mainly to the infrared final state singularities (compare fig. 16, fig. 17, fig. 19, 
fig. 21 and fig. 24), and moreover that although in some kinematic regions (high energy, low 
invariant mass and large rapidity) the NNLO correction to the hard cross-section may be 
important, the series converges sufficiently rapidly that NNNLO corrections are in practice 
usually small. This is reassuring. 

We also find that although in the resummed perturbation theory the number of in- 
frared logarithms increases by two for every extra power of a s , the effect of the extra 
logarithms is sufficiently benign that the hierarchy of the resummed perturbative expan- 
sion is not spoiled, and thus no further resummation is necessary. This again is due to 
asymptotic freedom, in the sense that it is only true when the coupling runs. 




(5.13) 
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A useful way to characterise this interplay between collinear and high energy resum- 
mation is to ask in which regions of the kinematic plane fig.l logarithms of Q 2 or logarithms 
of x are more important. If the former are dominant, we may expand the partonic cross- 
section H(N, Mi, M2) in powers of Mi and M2, keeping the full dependence on N, while 
if the small x logarithms dominate, we may expand in powers of N but keeping the full 
dependence on Mi and Mi- The relevant regions are thus characterised by the relative 
importance of factors of M and factors of N, or more specifically of the two integrals 

L,&t)= r^e* N ^^e^^)(M 1 + M 2 )L 2 (iV,M 1 ,M 2 ) = |L 2 (p,t), 
J_ ioo 2m J_ ioo 2m 2m at 

(6-1) 

These integrals may be computed for the gluon-gluon luminosity shown in fig. 12: the 
results are displayed in fig.l. Clearly when L' z ^> L z , high energy logarithms are relatively 
unimportant and the usual on-shell perturbative approximation to the hard cross-section 
is applicable. However when L' z C L z , the high energy logarithms are the most important, 
and we should use the off-shell perturbative expansion of the hard cross-section, perhaps 
with N = 0. It can be seen from fig.l that this is only true at very low scales, close to 
the initial boundary condition. More interesting is the intermediate region L' z ~ L z , in 
which resummation of both types of logarithm is necessary, and in which we may usefully 
expand both in N and in M. This region is important not only at high rapidity, but also 
in the central region when the invariant mass of the produced particles is not too high. 

The reason that the region dominated by logarithms of Q 2 is so much larger than that 
dominated by logarithms of x is very simple: when the coupling runs with Q 2 the variation 
of the luminosity with t = InQ 2 is much weaker than its variation with £ = In S/Q 2 , so 
factors of M are less significant than factors of N. This observation underpins the success 
of unresummed perturbation theory (which amounts to expansion in powers of M) in 
computing hard cross-sections in kinematic regions where one might naively have expected 
it to fail due to unresummed logarithms of x. 

We have shown a variety of estimates of the size of the high energy resummation 
effects compared to standard NLO perturbation theory. Our basic conclusion is that the 
effect of resummation in the partonic cross-section is an enhancement similar in size to 
that of a perturbative NNLO correction, while in the full cross-section the resummation 
of the gluon distribution produces an effect of similar magnitude but opposite sign. The 
net result is thus rather less than might be expected from NNLO considerations alone, 
since there is substantial cancellation. In fact we find that for a wide range of processes, 
each over a wide kinematic range, the net effect of resummation seems to be a suppression 
of between 5 and 10% (compare the blue curves in fig.ll, fig- 18, fig. 20 and fig. 25). This 
seems to suggest that the hard scale Q may not be the optimal factorisation scale for high 
energy processes. 

However it must also be remembered that the suppression of the gluon due to re- 
summed evolution is probably being overestimated in these calculations, since the initial 
distribution at 2GeV is kept fixed, rather than fitted to data (see for example the resummed 
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fits in refs. [26,38,39]). The true band of current uncertainty thus probably lies between 
these two extremes, i.e. between the solid and dashed blue curves in the figures. Thus 
for example in the total cross-section for hadronic S-production (fig. 20) the resummation 
corrections are as small as —5 ± 5% at the Tevatron, rising to 5 ± 10% at the LHC and 
20 ± 20% at a VLHC. These figures are still rather lower than suggested by the leading 
order calculations in ref.[14] and substantially smaller than the fixed coupling estimates 
in [15,16]. The effects on rapidity distributions may be rather larger, particularly at large 
rapidities: for B production at LHC resummation effects might be as large as ±15% at 
rapidities of 5 or so. Similar estimates should hold for inclusive jet cross- sect ions. 

For the gluonic contribution to Drell-Yan at LHC the effects are a little larger, since 
the infrared singularity is stronger. However for vector boson production the corrections 
are still modest since the hard scale is relatively large, though they are even so comparable 
to other sources of uncertainty. The predicted enhancement at large rapidities shown in 
fig. 24 is particularly striking. 

Refining these estimates into precise predictions is now straightforward: there are no 
longer any theoretical obstacles to computing cross-sections for hadroproduction processes 
correct to NLO which simultaneously resum all leading and next-to-leading logarithms of 
S and Q 2 . However there is still a lot of work to be done. 

Firstly it would be useful to have complete calculations of off-shell partonic cross- 
sections for a wider variety of hadronic processes: at the moment all we have are cross- 
sections for heavy quark production [16,27] and Higgs production in the m t — > oo limit 
[17]. The key partonic calculation for the future is clearly Drell-Yan and vector-boson 
production, both to confirm the conjectured structure of infrared singularities and provide 
a firm prediction for these benchmark processes. As explained above, it is important to 
match these cross-sections to the standard perturbative (on-shell) cross-sections, preferably 
by keeping the full N dependence when going off-shell. However for most purposes it 
is probably sufficient to treat the off-shellness perturbatively, keeping only the first few 
powers of M. The off-shell calculations are then complementary to the more usual fixed 
order (on-shell) calculations, offering a short cut to new results. 

Secondly, for more precise calculations we clearly need to include quark effects, par- 
ticularly in the high rapidity region where one of the partons is in the valence region. 
Including quarks in the resummed singlet evolution is no longer a problem [40], and in- 
cluding the contribution of initial state quarks to the resummed partonic cross-sections is 
also well understood [2,16,26]. 

Finally, for accurate resummed predictions it is necessary to produce resummed par- 
ton densities, fitted to data using resummed theoretical predictions. Previous experience 
[26,38] suggests that much of the effect of resummation in the evolution might then be ab- 
sorbed into the parton distributions, so without this ingredient resummation effects in the 
parton distribution functions are probably being overestimated. In order to obtain unbi- 
ased resummed parton distributions with sensible experimental uncertainty distributions, 
it will be necessary to use a statistical approach such as that currently being developed by 
the NNPDF collaboration [41]. 
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